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The formation of neutrino spectra in a supernova depends crucially on strength and inelasticity 
of weak interactions in hot nuclear matter. Neutrino interactions with nonrelativistic nucleons are 
mainly governed by the dynamical structure function for the nucleon spin density which describes 
its fluctuations. It has recently been shown that these fluctuations give rise to a new mode of energy 
transfer between neutrinos and nucleons which inside the neutrinosphere is of comparable or greater 
importance than ordinary recoil. We calculate numerically the spin density structure function 
in the limit of a dilute, non-degenerate medium from exact two-nucleon wave functions for some 
representative nuclear interaction potentials. We show that spectrum and magnitude of the energy 
transfer can deviate significantly from those based on the Born approximation. They are, however, 
rather insensitive to the particular nuclear potential as long as it reproduces experimental nucleon 
scattering phase shifts at energies up to a few tens of MeV. We also compare with calculations based 
on a one-pion exchange potential in Born approximation and briefly comment on their applicability 
near the center of a supernova core. Our study is relevant for numerical simulations of the neutrino 
spectra emerging from type-II supernovae. 

PACS numbers: 97.60.Bw, 13.15.+g, 14.60. Lm, 95.30.Cq 



I. INTRODUCTION 

The detection of roughly a dozen neutrinos from SN 
1987 A is in good qualitative agreement with the neutrino 
signal expected from the early cooling phase of a hot 
neutron star born in the center of the collapsed core of 
a massive star JjJ . It is therefore generally believed that 
type-II supernovae such as SN 1987 A are the optical 
counterparts of such catastrophic events. 

The formation of the spectra of neutrinos emitted from 
a type-II supernova takes place in a region where weak 
neutral current scattering and pair processes involving 
electron, /x, and r neutrinos and charged current creation 
and absorption of electron neutrinos on nucleons, nuclei 
and electrons cease to be efficient in keeping the neutri- 
nos in thermodynamical equilibrium with the medium. 
The interplay between (roughly) energy conserving scat- 
tering and energy changing reactions plays a crucial role 
in that respect |2|] . In previous studies of neutrino trans- 
port, the lowest order neutrino opacities in vacuum have 
been used. Neutral current scattering processes on nucle- 
ons and nuclei have been approximated to be elastic || . 
As a result, whereas the energy fluxes predicted for the 
three neutrino flavors turn out to be very similar , the 
effective temperatures are significantly higher for fi and r 
neutrinos compared to electron neutrinos which because 
of their more efficient energy exchange with the medium 
decouple from it further out. 

However, weak interaction rates in a medium differ sig- 



nificantly from those taking place in vacuum. On the 
one hand, the spin-dependent strong force between nu- 
cleons will establish spatial correlations of the density 
and the spin-density in the medium which can either en- 
hance or reduce average interaction rates. Many papers 
on weak interactions in neutron stars investigated these 
effects. However, they either applied the Landau theory 
of quasiparticles assuming a "cold" nuclear medium [ 4 £ ] 
or the authors focused on quasielastic scattering ||, 3 7 ] 
for which the energy transfer to the nucleons is smaller 
than the momentum transfer. 

On the other hand, a weak interaction transferring an 
energy u> to the medium is sensitive to the fluctuation 
power in density and spin-density at that frequency. For 
example, at finite density, the spin-dependent nucleon- 
nucleon interaction also causes the nucleon spins to fluc- 
tuate. This leads to a reduction of the average total axial- 
vector current neutrino scattering cross section compared 
to its vacuum value This effect is most important 

at the high temperatures pertaining in the first few sec- 
onds after formation of the hot neutron star. In addition, 
the nucleon spin fluctuations can, apart from recoil, im- 
ply an enhanced energy transfer between nucleons and 
neutrinos |h|. It is this effect which we are mostly con- 
cerned with in the present work because it could signifi- 
cantly change predictions of the neutrino spectra with a 
tendency to lower predicted effective temperatures of // 
and t neutrinos Q]. This is of some importance in view 
of new neutrino detectors such as Super-Kamiokande and 
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the Sudbury Neutrino Observatory, which have the capa- 
bility of measuring the neutrino spectra from nearby su- 
pernova events with much better statistics than is avail- 
able with the data from SN 1987 A. 

Near the surface of last scattering the neutrino opac- 
ities are governed mainly by the local nucleon spin den- 
sity. Within linear response theory, weak neutral-current 
interactions are then determined by the dynamical struc- 
ture function for the nucleon spin-density which describes 
its spatial and temporal correlations and is a function 
of energy and momentum transfer. For energy trans- 
fers that are larger than the typical spin fluctuation rate 
multiple nucleon-nucleon scattering is negligible, and, to 
lowest order in the nucleon-nucleon interaction, the spin- 
density structure function can be calculated from the ma- 
trix element for nucleon bremsstrahlung. This matrix 
element has been discussed in some detail in the litera- 
ture for a one-pion exchange (OPE) potential to 
lowest order in the pion-nucleon coupling, i.e. in Born 
approximation. 

However, the Born approximation is only applicable if 
at least one of the following conditions is fulfilled jl4| : 

\V\ « 

rriNa 

\V\ «-^, (1) 
rriMa 

where \ V\ ~ 100 MeV is the typical magnitude of the nu- 
clear interaction potential, a ~ 1 fm is its range, p is the 
nucleon momentum in the center of mass system, and 
tun is the free nucleon mass. The first condition is al- 
ways violated if the potential leads to bound states as 
for the proton-neutron interaction which gives rise to the 
deuteron bound state. The second condition translates 
into p > iriN and is therefore also violated for the non- 
relativistic nucleon momenta occurring in a supernova. 
We can therefore not expect that the Born approxima- 
tion is a reliable approximation to the dynamical nucleon 
spin-density structure function in a supernova. Neither 
is it obvious that any weak interaction rates calculated 
from it are reliable at the relatively low energies involved. 

The goal of this paper is therefore to compute the dy- 
namical nucleon spin-density structure function and re- 
sulting weak interaction rates beyond the Born approx- 
imation by using exact two- nucleon wave functions. To 
keep things numerically simple, we will restrict ourselves 
to spherically symmetric but spin-dependent two-nucleon 
potentials. Since a central potential conserves the total 
nucleon spin, the only contribution to inelastic weak pro- 
cesses (i.e. inelastic in the center of mass frame of the 
nucleons) will then arise from interactions of protons and 
neutrons due to their different weak coupling constants. 
We therefore have to deal with two nucleon species. Our 
approach takes into account in a consistent, unified way 
the free-free transitions 

ly + n+p^v + n+p (2) 



as well as the free-bound and bound-free processes in- 
volving the deuteron, 

i/ + d<-+i/+p + n. (3) 

The analogous processes involving neutrino pairs or ax- 
ions instead of neutrino scattering are described by the 
same dynamical structure function and can therefore also 
be treated by our formalism. 

The rest of the paper is organized as follows: In Sect. II 
we define the nucleon spin-density structure function in 
a form suitable for the case of proton-neutron interac- 
tions. The main formalism for computing this structure 
function from two-nucleon wave functions is presented in 
Sect. III. In Sect. IV the Born approximation is derived 
as a limiting case. Sect. V introduces a central potential 
which reproduces experimental data on proton-neutron 
scattering at energies below a few tens of MeV. The cor- 
responding spin-density structure function is numerically 
calculated for conditions around the neutrino sphere and 
compared with the Born approximation and calculations 
for an OPE potential. Finally, we summarize and con- 
clude in Sect. VI. 



II. THE NUCLEON SPIN-DENSITY 
STRUCTURE FUNCTION 

A. Definition and General Properties 

The interaction of interest here, neutrino-nucleon 
neutral-current scattering, is given by the Hamiltonian 

Hint = ~J= ^2 ~~ C A,i1b\ *l>i 

* i—n^p 

x^(l-7 S )^, (4) 

where Gf is the Fermi constant, ipi (i = n,p) and i/> v 
are the Dirac field operators for the neutrons, protons 
and neutrinos, and Cy,i and Ca,i are the dimension- 
less weak neutral-current vector and axial-vector cou- 
pling constants for protons and neutrons, respectively. 

Another possible type of weak process is the emission 
of axions [p.5fl . The corresponding interaction Hamilto- 
nian reads 

Hint = rj- ^2 Ca^il^lb^id^a , (5) 

a i—n,p 

where a is the axion field, f a the Peccei-Quinn scale, 
and the dimcnsionlcss coupling constants to neutrons and 
protons, C a ,i (i = n,p), are of order unity and depend 
on the specific axion model. 

In the limit of non-relativistic nucleons, only the axial- 
vector coupling contributes to inelastic weak processes. 
Within linear response theory these processes are then 
described exclusively by the dynamical structure func- 
tion for the nucleon spin density. In the following, we will 
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drop the index A or a in the nucleon coupling constants to 
neutrinos and axions, respectively, for notational simplic- 
ity whenever the result applies to both cases. To ensure 
a suitable normalization which will become clear below 
in Eq. (@), we define the structure function as the auto- 
correlation function of the weighted nucleon spin-density 



i—n,p 



(6) 



Here, (j>i(x) (i = n,p) is the non-relativistic field op- 
erator for protons and neutrons which is a Pauli two- 
spinor, r is the vector of Pauli matrices, and C 2 = 
C%Y n + CpY p is an average neutral-current axial weak 
coupling constant to the nucleons, weighted by the frac- 
tional neutron and proton abundances Y n and Y p . Defin- 
ing the Fourier transform in a normalization volume V as 
(T w (t,k) — V -1 / 2 f d 3 r e~ lk ' r cr w (t,r), one can then de- 
fine the nuceon spin-density structure function |lj|,[l^,[l6| : 

4 r+°° 

S„(u;,k) = — dte wt {a w (t,-k)-a w (0,-k)). (7) 
OUb J _ OQ 

Here, (u>, k) is the four-momentum transfer to the 
medium, n b is the baryon number density and the ex- 
pectation value (• • •) is taken over a thermal ensemble at 
the medium temperature T of medium states normalized 
to unity. 

Relativistic neutrinos and possibly axions will have 
typical energies of order 3T but are in general not in 
chemical equilibrium with the medium. Weak interac- 
tions such as neutral current neutrino scattering and pair 
processes and axion emission thus probe the spin-density 
function typically at thermal energy-momentum trans- 
fers. Since the momenta involved in the nucleon-nucleon 
interactions are much larger than the thermal momenta 
of relativistic particles, we will often employ the long 
wavelength limit, S a (u>) = limk— >o S a {u}, k) for which we 
define the dimensionless quantity S a (x) = TS a (xT). In 
this limit, integration of Eq. (|^) over u> yields the sum 
rule 



+oo 



did 

2^ 



S a (u) -1 = N„ 



3n b V 



(8) 



where we wrote the spatial integral J d 3 r cr w (t, r) = 
J2i a 'i,w Here ) °»,u> = °"« diag(C p , C n )/C, where er; are 
the spin operators of the individual nucleons, and the ma- 
trix diag(C p , C n ) acts in isospin space. In Eq. (||) N a de- 
scribes correlations among different nucleon spins. In the 
absence of such correlations f*™(dw/2ir)S a (u)) reduces 
to 1 which motivated the introduction of the weighted 
spin operator Eq. (||). 

We formally introduce the complete set of eigenfunc- 
tions \n) of the total Hamiltonian H of the nuclear 
medium, H \n) = Lo n \n), where Lo n are the corresponding 



energy eigenvalues. By inserting the identity operator 
/ = \n) (n\, between the spin operators, Eq. (Q) can be 
rewritten into 



&7T 1 

S^c^k) = - — — y e~ Un/T 5(u) + u) n -u> m ) 
drib Z ' 

x |(n|,r„(0,k)|TO)| 2 , 



(9) 



where Z = ^ n e~^ n ^ T is the partition function. This 
form will be useful later and it shows that the structure 
function satisfies detailed balance, 



5 CT (w,k)=5 a (-w,-k)e w / T . 



(10) 
for 



It is therefore sufficient to know the function, e.g. 
positive energy transfer to the medium, lu > 0. 

Up to now no approximations have been made with re- 
gard to the nucleon-nucleon interactions which determine 
the nonperturbative though unknown structure function 
5 CT (cj,k). From now on we will make the assumption 
that only two-nucleon forces are present. The Hamilto- 
nian then has the form 



2m N 



131 ' 



ii 17 j) 



(11) 



where is the radius vector between nucleon i and j, Pi 
is the nucleon momentum, V(ry, &i, <Tj) is the spin de- 
pendent two-nucleon interaction potential, and the sums 
run over all nucleons. The most general two-nucleon po- 
tential in the non-relativistic limit can be written as [T^ ] 

V(t, 0-1,0-2) = U(r) + U a (r)<Ti • <x 2 + U T (r)T 12 (12) 
+P T [U T {r) + UHr)*i ■ (7 2 + U^(r)T 12 ] 

where r = ri2, r = |r|, f = r/r, P T is the isospin exchange 
operator, and the tensor operator is given by 



T12 = 3 <T\ f CT 2 f — (J\ ■ CT2 



(13) 



Useful information about structure functions is con- 
tained in their moments of which Eq.(|^) is an example 
for the lowest one. The next higher moment is given by 
the so called f-sum rule which is often discussed in the 
literature in the context of the density structure function 
and for spin-conserving interactions |ll|. In Ref. |l6| ] 
we derived a generalized f-sum rule for the spin-density 
structure function for one species of nucleons interact- 
ing via spin-dependent forces of the form Eq. ( |l2|) in a 
non-degenerate medium: 



— uS a (uj) = — 

00 2?r n b V 



H 7 



(14) 



Here, Ht is the part of the total Hamiltonian involving 
the tensor operator Ty . In the present paper we will 
consider both neutrons and protons but assume a central 
two-nucleon potential, i.e. absence of tensor contribu- 
tions. The f-sum rule is then modified to 
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3n b V 



c 



TO, (15) 



where H™ p is the spin-dependent central part of the to- 
tal Hamiltonian which contributes to neutron-proton in- 
teractions. Note from Eqs. and ( |l4| ) that for only 
one nucleon species a tensor interaction is required to 
give a non-trivial spin-density structure function. This 
is because the central part of the interaction conserves 
the total nucleon spin and thus does not contribute to 
its fluctuations. In contrast, for two nucleon species, a 
central spin-dependent proton-neutron interaction is suf- 
ficient for a non-trivial structure function as long as the 
neutral-current axial weak coupling constants for protons 
and neutrons are different [see Eq. dll])]. We stress, how- 
ever, that the actual (positive) value of the f sum depends 
on all interaction terms via the states entering the ther- 
mal average. 



B. Relevance for Weak Interactions 

The differential axial-vector-current neutrino-nucleon 
cross section is determined by the dynamical nucleon 
spin-density structure function S (T {u,'k), taken at the 
difference of initial and final neutrino four-momentum 
(ei,ki) and (e 2 ,k 2 ) via 



do A = GrC 2 



3_cosff d 3 k 2 



(2tt) 3 



(16) 



where 6 is the angle between ki and k 2 . In our con- 
vention, the neutral current axial-vector contribution to 
neutrino scattering rates on the ensemble of all nucleons 
is ribdcFA- 

The axion emission rate per volume, Q a , is governed by 
a structure function S aya which is obtained from Eqs. (|^), 
(0) by substituting C t — * C a>i (i = n,p), 



C 2 a n b 
i^YPa Jo 



(17) 



where C 2 — C% Y n + C% p Y p . We have assumed an 
isotropic medium such that S^cjjk) only depends on 

fc= ' k l- 

Various quantities relevant for neutrino diffusion are 
determined by the spin-density structure function. For 
the remainder of this section, we assume a Maxwell- 
Boltzmann distribution at temperature T v for the neu- 
trinos. Furthermore, we make use of the normalization 
given by Eq. (||). Contributions form spin-spin correla- 
tions represented by N a are mainly induced by the pres- 
ence of nucleon bound states and by the Pauli exclu- 
sion principle which becomes important in a degenerate 
medium ]l!j. Both effects are small in the post-collapse 
phase of a supernova in which we are interested. We 
therefore assume N a <C 1 in Eq. (|J). The average en- 
ergy transfer per collision in a dilute medium can then 
be written as |10L|9|1 



(Ae) 
T 







°° dx ~ ( J3x 2 
-S„{x)lx+ — 



~L2~ 



(18) 



with j3 = T/T V . This should be compared to the average 
energy transfer by nucleon recoils KG], 



(As) 



recoil 



30(/3- 1) T 2 



fl 2 



m N 



(19) 



Another interesting quantity is the reduction of the av- 
erage total axial-vector current scattering cross section 
(a a) [see Eq. (pi|)] in the nuclear medium First we 

note that a term of the form AS(u>) in S a (u>) corresponds 
to a total elastic scattering cross section 



o-ez(ei) 



3A 
8^ 



G'2. r~i2 2 



(20) 



For rib —> Q there are no spin fluctuations and corre- 
lations, and Eq. (@) implies S a (ui) — 2ttS(uj) and thus 
(T = (9 / n)C A G 2 F T 2 for the thermally averaged cross sec- 
tion. In Ref. [§ we obtained the expression 



6(<t a ) 



cto 



= N„- 



(21) 



ri*<.> 

Jo 27r 



l + x + — 



which again holds in the dilute medium. The physical 
quantities discussed here will be calculated for the su- 
pernova environment in Sect. V. 



III. BEYOND THE BORN APPROXIMATION 

A. Classical versus General Quantum Result 

In the limit |u;| <C T the nucleon spin can be treated as 
a classical spin s being changed abruptly by some random 
amount As in a typical nucleon-nucleon collision event 
which takes place on a time scale ~ 1 /T and thus appears 
to be "hard" . In this case we expect [ElJ 



co 2 + r 2 /4 ' 



(22) 



where the spin fluctuation rate r CT is related to the colli- 
sion rate r co ii by 



v„ = 



((As) 2 ) 



coll • 



(23) 



Note that the spin fluctuation rate suppresses the uj~ 2 
bremsstrahlung spectrum which otherwise would vio- 
late the existence of the normalization Eq. (||). This 
is known as the Landau-Pomeranchuk-Migdal (LPM) ef- 
fect MM. In previous work pl^Mpf it has been 



4 



discussed how the Lorentzian shape Eq. ( |22| ) might in- 
fluence weak interaction rates at high densities where 
r CT > T. The high density behavior of the spin-density 
structure function can also influence limits on the axion 
mass pt| . 

For \u>\ ^> T a multiple scattering effects can be ig- 
nored and the spin-density structure function can be 
computed by a quantum mechanical treatment of two- 
nucleon scattering. From the generic tu~ 2 divergence of 
all bremsstrahlung processes for u> — > one expects the 
general form || 



%^/T)x(^ /T fo ^<°' (24) 



I 



for uj > 0, 



where s(x) is a nonsingular even function with s(0) = 1. 
The specific shape of s(x) for x > 1 depends on the 
nucleon-nucleon interaction potential Eq. (|l2|), and its 
calculation for realistic interaction potentials is the main 
goal of this paper. Comparing Eqs. (|2^) and ( |24| ) in 
their common range of validity, r CT <C <C T, shows 
that the coefficient T a of the bremsstrahlung divergence 
in Eq. (^) can be interpreted as a nucleon spin fluctu- 
ation rate and that the classical limit of hard collisions 
corresponds to s(x) = 1. The existence of the f sum 
Eqs. (fhf), ( |l5| ) shows that s(x) has to decrease for large 
x due to quantum corrections. 



B. Exact Treatment in the Limit of High Energy 
Transfers 

For lo > where scattering involving more than two 
nucleons is negligible, we can numerically compute two- 
nucleon wave functions from a given nucleon-nucleon in- 
teraction potential and use it in the general expression 
Eq. (^). For a central potential the eigenfunctions for 
the relative motion in the proton-neutron center of mass 
system 



IP) 



\p,l,m,S) = R pls (r)Y lm (n)\S) 



(25) 



are characterized by the quantum numbers for the ra- 
dial momentum, p, the orbital angular momentum, I 
and m, and the total spin, S, where R p is(r) is the ra- 
dial wave function and Yi m (fl) are the spherical harmon- 
ics. The corresponding energy eigenvalues up have a 
(21 + 1)(2S + 1) fold degeneracy. The Pauli exclusion 
principle then also determines the isospin to I = | + 
(— 1)' — S 1 ). Assuming an isotropic medium and us- 



ing <t w (0, k) - V" 1 / 2 ^^ 



diag(C p , C n )/C with 



a normalization volume V = 1/rib, after some algebraic 
manipulations we obtain 



(26) 



/T-rn N [uj+u P -uj Q -k 2 /{Am N )\ 2 /{Tk 2 ) 



P 



Cn 



where Zcm — J2p e UJp / T . For k — > this expression 
transforms into 



(27) 



x E E( 2Z + x ) E ^ P/T <5(^ + - wq) 

P,g I S=0 

2 

2 D* 



X 







drr RUs( r ) R qi(i-s)(r) 



where we have made use of the orthogonality of the sys- 
tem of eigenfunctions which are supposed to be normal- 
ized to unity. In practice one constructs bound and scat- 
tering states of the stationary radial Schrodinger equa- 
tion within a finite spherical volume of radius r max and 
computes the matrix elements appearing in Eq. (p7|). For 
a nucleon interaction potential that is not radially sym- 
metric, the eigenstates \P) = \pJV) are characterized 
by the total angular momentum J and parity V and are 
superpositions of orbital angular momentum eigenstates. 
With this modification, Eq. (|2^) still holds but we will 
not pursue this more complicated case here which would 
lead to coupled radial equations for the corresponding 
radial functions R p j-p. 

Since we neglect interactions among more than two 
nucleons, our formalism does only account for neutrons, 
protons and deuterons. Higher nuclei such as helium are 
not included. In this sense, strictly, rib, Y pi and Y n have 
to be interpreted within the ensemble of neutrons, pro- 
tons and deuterons only. Nuclear statistical equilibrium 
shows that in practice this does not make a big difference 
in our situation where Y p <C Y n . Keeping this in mind 
we can now calculate Zcm analytically. Around the neu- 
trinosphere the nucleons are at best mildly degenerate. 
We therefore assume a Maxwell-Boltzmann distribution 
f(p) = e~ p /( 2 ^ T ) for the unbound proton-neutron states. 
Here, p — s 1 / 2 /2 is the nucleon momentum in the center 
of mass system, expressed in terms of the squared cen- 
ter of mass energy s (excluding the nucleon rest mass), 
and \i = mjv/2 is the reduced nucleon mass. Taking into 
account the spin degrees of freedom we then have 



Zcm = 3e 



n b \2tt ) 



3/2 



(28) 



where ~ 2.2 MeV is the deuteron binding energy (the 
deuteron has 5=1). The degree of dissociation, i.e. the 
fractional abundance of unbound states is then 



fu 



3/2 



-e d /T 



(29) 



As a consequence, S a (u)) from Eqs. ( |26| ) and ( J27| ) does 
not exhibit a simple linear scaling with the nucleon den- 
sity rib, except for the dilute limit, rib —> 0, fa —> f. It 
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can be seen that the numerator in Eqs. ( |2q ) and ( p7| ) is 
independent of n b and the density dependence of S a thus 
exclusively stems from Zqm- 

In the limit of zero temperature, /„ — > and only 
the deuteron bound state will be populated in thermal 
equilibrium. In this limit, Eq. ( jig ) describes the cross 
section for the weak neutral current deuteron break up 
process, Eq. (||). Integration over the phase space for the 
outgoing neutrino yields the total cross section 



3G|i (Cp — C n ) z 



16tt 2 



Y 

doj(ei — u) 2 lim S a (ui) 

IT~>0 



(30) 



for incident neutrino energy e\. 

Another instructive limiting case is the absence of spin 
flip interactions. Scattering on protons and neutrons 
then has to be added incoherently with the states \P) 
now being plane waves. Eq. (|2(f ) then reduces to the 
ordinary recoil expression 



S, 



1 /2 



(31) 



Let us now get back to the general case. In agreement 
with Eq. (|l5|) and Ref. [Q only proton-neutron scat- 
tering contributes to Eq. Q27[), and only if the neutral- 
current axial weak coupling constants for protons and 
neutrons are different. Since the total spin is conserved 
by a central potential, the spin-density structure func- 
tion is governed by the fluctuations of the difference of 
the proton and neutron spin, <j p — cr n . To compare with 
the general results Eqs. ( p2| ) and (|24|), we study the cor- 
responding spin-flip cross section which is defined as 



((<x p - tr n ) ) 



(32) 



Here, a np (s) = ^ ; <7 rap (/,s), where the average total 
proton-neutron scattering cross section in angular mo- 
mentum state / is given by 

An 

(T np (l,s) = —(21 + 1) [3sin 2 (5 ;!l ( S )+sin 2 ^o(s)] (33) 

in terms of the phase shifts Sis(s). The latter are defined 
by the asymptotic behavior 

R p is(r) oc sin [(2^)^^ - Ztt/2 + 5 is (8(ilu p )\ , (34) 

of the scattering states uip > for (2/iujp) 1 ^ 2 r 3> I- The 
nucleon spin flip rate is now just defined as 



T s / = Y p Y„n b 



J +c °dpp 2 f(p)(p/^a sf (4p 2 
jo + °° dpp 2 f(p) 



(35) 



where pj \x is the relative velocity of proton and neutron. 

As can be seen from phase shift analysis, the spin flip 
cross section Eq. ([}2]) is 

<V(«) = IT + X ) M n ^,i( s ) - sm<5 ; , (s)] 2 . (36) 

OS 



Note that this vanishes if the phase shifts for S = 
and S = 1 are equal, as expected. Given a s f, one can 
compute T s f from Eq. (^5|) and compare it with the Born 
approximation to be discussed below and with Eq. J2~i|). 
This will be done in the following two sections. 



IV. THE BORN APPROXIMATION 

By expanding the unbound states \P) into plane waves 
within first order perturbation theory and inserting the 
result into Eq. (^6|), one obtains the spin-density struc- 
ture function in Born approximation. For u> > 0, the 
result in the long wavelength limit is 



^ orn H 



^YYn b 



J +co dppf( P )j^;;dkk\u^(k)\ 2 

p + OC 



(37) 



Jo °° dpp 2 f(p) 

where fc m ax,min = (p 2 + 2/xw) 1 / 2 ±p, and U™ p (k) is the 
Fourier transform of the coefficient of cr p ■ <j n in the 
proton-neutron interaction (er p and er„ being the proton 
and neutron spins). Only the relative motion between 
proton and neutron influences Eq. (|37]) because neither 
energy nor momentum can be transferred to the center 
of mass motion in the long wavelength limit. In contrast 
to Eq. (§7|), S% ora (u>) scales linearly with n b . 

In Ref. ft we considered one species of nucleons cou- 
pling to a classical, external scattering center via an in- 
teraction of the form Eq. (|l^) where one of the spins 
was replaced by a classical spin s associated with the ex- 
ternal scatterer. The result for the spin-density structure 
function in the long wavelength limit in Born approxima- 
tion is very similar to Eq. (J37|) for the case of a central 
two-nucleon potential and a medium of protons and neu- 
trons with different neutral-current axial weak coupling 
constants. 

In Born approximation, the spin flip cross section 
Eq. (p2l) evaluates to 



20 a 2 f sl ' 2 



(38) 



where ((er p — cr n ) 2 ) = 3/2 was used. Comparing 
Eqs. (Eg), ©, and pi) yields 



27 / Cp C n 
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c 



pBorn 



(39) 



G 



i.e. the spin fluctuation rate in S a (u)) and the average 
spin flip rate indeed agree within a factor of order unity, 
apart from the factor [(C p — C„)/C] 2 involving the weak 
coupling constants which results from our specific defini- 
tion of S a . 

As an example, we consider the usually adopted OPE 
potential which is a good approximation to the nucleon- 
nucleon interaction for distances greater than the inverse 
pion mass m % . With / ~ 1 being the pion-nucleon cou- 
pling constant, its Fourier transform is 



2/V ( gl • k) (Q-2 ■ k) 

iL2 _|_ 7^-1 2 



VbpE(k,<T 1 ,er 2 ) = - I 



(40) 



More generally, as can be seen from Eq. 



oc to 



~ 3 / 2 ~ r for w 



OC 



if \U a {k)\ 



]), one has 

cx k~ r for 



5 < 

/c — ► oo , corresponding to existence and square integrabil- 
ity of the (r — 2)th derivative of the interaction potential. 
This should hold independently of the Born approxima- 
tion which is viable in the limit of high energies as we 
will see now. 



V. A NUMERICAL MODEL FOR THE 
SUPERNOVA ENVIRONMENT 

We first note that S a {u) from Eqs. © and © is 
proportional to the dimensionless factor 



and it clearly has a tensor contribution. The spin- 
density structure function corresponding to this potential 
has been calculated in Born approximation Jll|. Trans- 
lated into our notation, the contribution from proton- 
neutron scattering takes the form of Eq. ( gjj ) with s(x) = 
s(x)/s(0) given by the function 



s(x) 



dv[v{v + x)] 1/2 e- v 



(41) 



(5C 2 + 3C 2 ) si(v, x) + 2 (C 2 + C 2 _) s 2 (v, x) 



(6C 2 + + 2C 2 _)s 3 (v,x) 



where C± — (C p ± C ra )/2 and 



dz < 



2v+x-2z[v(v+x)] 1/2 
2v+x-2z[v(v+x)] 1 / 2 +y 

(2v+x) 2 -iv(v+x)z 2 
{2v J r x J ry) 2 —A:v{v J rx)z 2 



„ (2v+x+y) 2 -±v(v+x)z 2 



i = 1, 
i = 2, 

i = 3, 



(42) 



with y = rn^^misiT). Furthermore, the contribution to 
the nucleon spin fluctuation rate T a is 



r 



= -Y Y 

1 r) J r, 



5(0), 



o\OPE — g 1 p 1 n ~Q2^ A 



(43) 



with 



2 n b 



T l/2 



5/2 



= 8.6 MeV pi 3 Tj 



1/2 



Here, a* = (/2m A r/m 7r ) 2 /47r w 15, p i3 
with p the density, and T w = T/lOMeV 



(O /10 13 gcm 3 



Note that 



S^qp E (uj) oc cj 3 / 2 for a; — > oo and thus violates f sum 
integrability. As was explained in Ref. jl6| , this is caused 
by the unphysical behavior of the OPE potential for 
r — > that leads to an |VbpE(k)l which for k — > oo is 
asymptotically constant [see Eq. (|40|)]. 



J — 1 p 1 n 



c 



(Cp C n ) 



y y 

± p 1 n 



C 2 Yp 



C 2 Y n 



(45) 



which describes its compositional and coupling constant 
dependence for fixed rib and T. For processes involving 
only protons or neutrons this factor would be replaced 
by Y 2 and Y 2 , respectively. Since interaction rates are 
proportional to nbS a by definition, Y is a rough mea- 
sure of the contribution of proton-neutron scattering to 
weak neutral current inelastic interaction rates. For the 
neutrino-nucleon coupling in a nuclear medium we will 
adopt C A . P ^ 1-09, and C A ,n ^ -0.91 (ij), so that 
Y ~ 0.5 for Y p ~ 0.1. 

In the following we are interested in the environment 
given inside but not far from the neutrinosphere in a su- 
pernova. For the rest of this section we choose the repre- 
sentative numbers p = 10 13 gem -3 , T — 8 MeV, and T„ = 
10 MeV. For these parameters, /„ ~ 0.36, corresponding 
to a fractional deuterium abundance Yd — (1— fujYp. Nu- 
clear statistical equilibrium involving higher nuclei gives 
values that are within 20 - 30% of this if Y p < 0.2. We 
stress again that due to the presence of bound states 
the spin-density structure function calculated by phase 
shift analysis and weak interaction rates computed from 
it do not exhibit a simple scaling behavior with den- 
sity and/or temperature, as discussed below Eq. (gg). 
Our choice represents a typical case of interest for the 
neutrino-nucleon energy transfer. 

For the proton-neutron interaction potential Vg P (r) 
for total spin S we chose the following Gaussian poten- 
tials [such that U^ p = V" p - V£ p and U np = {V™ p + 



(44) Vz P )/2 in the notation of Eq. 



V np (r) = -33.6e- (r/1 - 77fm)2 MeV 
y"P( r ) = _84.7 e -(»-/i-36fin) a M eV. 



(46) 



Its strengths and ranges were fit to reproduce the exper- 
imental values for the scattering lengths as and effective 
ranges which determine the low-energy expansion 

of the phase shifts 5q,s @ : 
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FIG. 1. The s-wave proton-neutron scattering cross sec- 
tion predicted by the potential Eq. ([ifj]) (thin solid line), and 
measured (thin dotted line) as a function of the laboratory 
kinetic energy Ei a b = s/2mjv. Also shown for this potential 
is the spin flip cross section in Born approximation [Eq. (38), 
thick dashed line] and from phase shift analysis [Eq. (36), 
thick solid line]. The inverted resonance in the latter curve 
at £i a b/2 ~ 2.2 MeV stems from the deuteron bound state. 



1 v 

cot 6 ,s(s = 4p 2 ) = + £ r cS 3 . (47) 

pas 2 

As a result, the s-wave proton-neutron scattering cross 
section predicted by Eq. fliq ) agrees with the experimen- 
tal one to within less than 5% in the laboratory energy 
range between and ~ 20 MeV (see Fig. [l]). In addition, 
the energy of the bound state resulting for S — 1 co- 
incides with the deuteron binding energy within 5%. A 
central potential describes the deuteron rather well since 
the contribution of the D state to the bound state wave 
function is only about 6%. Finally, we have compared the 
numbers for the weak neutral current deuteron break up 
cross section resulting from Eq. ( |3C| ) with calculations in 
the literature p7| . In the energy range between ~ 5 MeV 
and 40 MeV we found agreement to within about 10%. 
This serves as a further check for the correct normaliza- 
tion of our calculation. 

Also shown in Fig. [l] is the spin flip cross section as cal- 
culated from the potential Eq. ( fl6|) both in Born approx- 
imation [Eq. (p8|) ] and numerically from the phase shifts 
[Eq. (p6|)1. It is clearly seen that the Born approxima- 
tion is far from being good. The integrated spin flip rate 
Eq. @, r s/ ~ 0A6Y p Y n MeV, differs by about a factor 
2.5 from the Born approximation r^ orn ~ 0.21Y p Y n MeV. 
This is also reflected by the nucleoli spin-density struc- 
ture function calculated from Eqs. (37) and ( p7[ ) for 
k —>■ 0, as shown in Fig. g. The phase shift analysis was 



FIG. 2. The contribution to the dynamical nucleon 
spin-density structure function S a (ui)/Y from proton-neutron 
scattering as a function of u) in the long wavelength limit 
k — > 0. Shown is the Born approximation [Eq. (pj), thick 
dashed line] and the result from phase shift analysis [Eq. (^) , 
thick solid line] for the proton-neutron interaction potential 
Eq. @ and the estimate Eq. @ with Eqs. @-0) for 
an OPE potential in Born approximation (thin dashed line). 
The resonance at uj ~ 2.2 MeV from the deuteron binding 
energy is clearly visible in the thick solid line. The small 
wiggles on this curve are caused by the finite numerical reso- 
lution of the energy eigenvalues of the scattering states. For 
uj < F a ~ 1.5 MeV, multiple scattering effects start to become 
important and tend to suppress S a (u)/Y below the values 
shown here. 



performed by computing the radial eigenfunctions up to 
some maximal orbital angular momentum Z max = 3 above 
which they are close enough to the free eigenfunctions to 
make a negligible contribution to Eq. (^7|). To achieve a 
sufficient resolution in the energy range of interest, about 
500 eigenfunctions had to be computed. We verified that 
the resulting S a (oj) satisfies the f sum rule Eq. @ to 
within 10%. Note, furthermore, that the Born approxi- 
mation and the phase shift calculation of the quantities 
shown in Figs. |l], and ^| converge at high energies where 
the second condition in Eq. ([y) is asymptotically satis- 
fied. We also verified that for a weak interaction potential 
satisfying the first condition in Eq. (Q), the Born approx- 
imation agrees well with the phase shift analysis over the 
whole energy range as expected. 

For comparison, Fig. || also shows Eq. (|24|) with the 
perturbative expressions Eqs. (41 -|44|) for the proton- 
neutron scattering contribution to S a based on the OPE 
potential (thin dashed curve). After all, this curve re- 
produces the general normalization of the spin-density 
structure function quite well, but it cannot reproduce the 
quite prominent deuteron resonance. The non-vanishing 



8 



pion mass is taken into account in this curve and sup- 
presses it by roughly a factor 2 compared to calculations 
neglecting the pion mass. Note the steepening at high 
iv of the curves for the potential Eq. ( flrf ) in contrast to 
5^ope( w ) which guarantees or violates f sum integrabil- 
ity, respectively. 

The quantities of interest for neutrino diffusion as dis- 
cussed in Sect. II are shown in Table |[ In calculating the 
cross section reduction Eq. (|2l| ) we neglected the term N a 
involving spin-spin correlations. It is easy to see that the 
5=1 deuteron bound state yields the contribution 
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N a = 



C A,pC A,n Yd 



a 2 



(48) 



which implies a further reduction of (a a) because of the 
opposite sign of Ca, p and Ca,u- This corresponds to 
the fact that the cross section for elastic scattering on 
deuterons is significantly smaller than that on free nu- 
cleons. Also note that Eqs. (18) and ( |l9| ) imply that 
(Ae) oc T v - T for \T V - T\ <C T. Finally, we observe 
that T CT /r s/ ~ T^ oln /T^ ln , i.e. the relation Eq. @ 
between the spin fluctuation rate T„ appearing in S a (u>) 
and the spin flip rate T s f also holds beyond the Born 
approximation. 

IV Y [MeV] (Ae) /(T„ - T) S (a A ) /(a Y) 



recoil 

OPEJ3orn 
Eq. (Get), Born 
Eq. (46), exact 



3.6 
0.66 
1.5 



0.32 
0.18Y 
0.029F 
0.15Y 



-0.14 
-0.018 
-0.076 



TABLE I. Comparison of the quantities shown in the 
top row for the cases specified in the first column for pro- 
ton-neutron interactions. The first entry in the first column 
is the proton-neutron potential used, and the second entry in- 
dicates whether the Born approximation ["Born", according 
to Eq. ( pi7| )] or the exact two-nucleon wave functions ["exact", 
according to Eq. (p?j)] have been used. For the OPE case, 
C p = 1.09 and C n = —0.91 have been assumed. Spin-spin 
correlations have been neglected in S {a a) /(To (see text). 



Clearly, the Born approximation for the potential 
Eq. (^6|) does not give a good estimate to any of these 
quantities. However, the numbers based on the OPE po- 
tential in Born approximation give a reasonable estimate 
to most of the integrated quantities in Table Q. In par- 
ticular, these results confirm that the average inelastic 
neutrino-nucleon energy transfer (Ae) is indeed compa- 
rable to the recoil energy (Ae) reC oii, as suggested by cal- 
culations employing the Born approximation for the OPE 
potential |g,|. This energy transfer is, however, differ- 
ently distributed with a much longer tail to high energy 
transfers, as can be seen in Fig. ||. 

We have furthermore checked that the results for So- 
calculated from the phase shift analysis Eq. ( p7| ) are in- 
sensitive to the detailed shape of U™ p (r) as long as it 
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3. The energy 

transfer distribution S<?(xT) [l + f3x/2 + (fix) 2 /12] , resulting 
from proton-neutron scattering, i.e. the integrand of Eq. jl^ ) 
divided by ui, as a function of x — u)/T. Here, the thick solid 
and dashed lines were obtained from the phase shift analysis, 
Eq. (pTj), and from the Born approximation, Eq. j37j), respec- 
tively, for Y = 1 and k — > 0. For comparison, the thin solid 
line line was obtained from the recoil energy transfer, Eq. (|jl|), 
for a typical thermal momentum transfer k = 37^, and the 
thin dashed line is from Eq. (^) with Eqs. (0-0) for an 
OPE potential in Born approximation. For x < Y a /T ~ 0.2, 
multiple scattering effects start to become important and tend 
to suppress the energy transfer distribution below the values 
shown here. 



reproduces the experimental phase shifts in the corre- 
sponding energy range. In particular, properties of the 
potential at short distances r influence So-(u>) only for 
u> = p 2 /win > l/(ni]yr 2 ). For the conditions near the 
neutrinosphere it is therefore sufficient that the potential 
reproduces nucleon-nucleon scattering up to a few tens 
of MeV. 

Towards the center of the hot neutron star, at densities 
around nuclear density and T ~ 30 — 50 MeV, predictions 
for the quantities shown in Table | by the OPE potential 
in Born approximation are about 10 times higher than 
corresponding predictions based on the potential Eq. d46| ) 
for which Born approximation and phase shift analysis 
become rather similar. This shows that in this environ- 
ment weak interaction rates become quite sensitive to the 
short distance behavior of the two-nucleon interaction 
potential which is different for these two potentials. This 
can have important ramifications for neutrino opacities 
and axion emissivities in the supernova core that are usu- 
ally based on these OPE calculations [ pd|p2| . Whereas 
calculations assuming an OPE potential should be a rea- 
sonable approximation in the context of a "cold" neu- 
tron star this is not necessarily the case for the much 







higher thermal energies involved in a hot protoneutron 
star. Neutrino opacities govern the cooling time scale of 
the protoneutron star pj| while axion emissivities deter- 
mine axion mass bounds based on supernovae [^6| . Apart 
from taking into account many-body effects such as mul- 
tiple scattering |24|,[l(J[l6| a more reliable calculation of 
these quantities thus requires to use nuclear potentials 
that fit nucleon-nucleon scattering data also at energies 
above a few tens of MeV to ensure the correct small dis- 
tance behavior. We leave that to a separate study. 



VI. SUMMARY AND CONCLUSIONS 

We have discussed weak axial-vector current interac- 
tions involving nucleons in hot non-degenerate nuclear 
matter at temperatures around 10 MeV and densities of 
a few percent of nuclear saturation density, i.e. for con- 
ditions given in the vicinity of the neutrinosphere in a 
type-II supernova. To describe such interactions in the 
limit of non-relativistic nucleons we adopted the struc- 
ture function formalism for the nucleon spin-density. Our 
special emphasize was on the energy transfer between the 
weak "probe" and the nuclear medium which is induced 
by the nucleon spin fluctuations caused by the spin- 
dependent nucleon-nucleon interactions. To lowest order 
in the nucleon-nucleon interactions, i.e. in Born approx- 
imation, this is represented by nucleon brcmsstrahlung. 
We have shown, however, that the Born approximation 
is in general not a reliable estimate for these effects. As 
an alternative, we have performed computations using 
exact two-nucleon wave functions for a spherically sym- 
metric two-nucleon interaction potential that was fit to 
experimental data. In this case, only proton-neutron in- 
teractions contribute to inelastic weak neutral-current in- 
teractions with nucleons. We compared our calculations 
with results for the corresponding contribution based on 
the usually adopted OPE potential in Born approxima- 
tion. One gets rather good agreement for most integrated 
quantities such as the reduction of the total neutrino scat- 
tering cross section and the average energy transfer in a 
scattering event. In particular, we confirm that the lat- 
ter is comparable to the recoil energy, as suggested by 
the results for the OPE potential in Born approximation. 
In contrast, differential quantities such as the distribu- 
tion of energy transfers deviate significantly from predic- 
tions based on the Born approximation. Our calculations 
specifically predict that the energy transfer peaks at the 
deuteron binding energy. Resonances from higher nuclei 
which have not been taken into account here such as he- 
lium are probably less important for T < 10 MeV because 
they would appear at higher energies and contribute less 
to the energy transfer distribution (see Fig. ||). 

The formalism presented here can be extended to two- 
nucleon potentials that are not spherically symmetric and 
to finite momentum transfer [see Eq. (p6[)]. Our results 
might have a significant impact on the formation of neu- 



trino spectra from type-II supernovae. A quantitative 
understanding will, however, require detailed numerical 
simulations. Finally, we demonstrated that weak interac- 
tion rates in the hot supernova core are sensitive to the 
small distance behavior of the nucleon-nucleon interac- 
tion potential which is not well described by the usually 
adopted OPE potential. This should be taken into ac- 
count in future investigations. 
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